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Abstract: Non-uniform black strings (NUBS) are studied by the large D effective theory 
approach. By solving the near-horizon geometry in the 1/D expansion, we obtain the 
effective equation for the deformed horizon up to the next-to-next-to-leading order (NNLO) 
in 1/D. We also solve the far-zone geometry by the Newtonian approximation. Matching 
the near and far zones, the thermodynamic variables are computed in the 1/D expansion. 
As the result, the large D analysis gives a critical dimension D* ~ 13.5 at which the 
translation-symmetry-breaking phase transition changes between first and second order. 
This value of D* agrees perfectly, within the precision of the 1 /D expansion, with the result 
previously obtained by E. Sorkin through the numerical resolution. We also compare our 
NNLO results for the thermodynamics of NUBS to earlier numerical calculations, and find 
good agreement within the expected precision. 
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1 Introduction 


Black holes in the dimension D > 4 have totally different aspects than in D = 4. A remark¬ 
able property of higher dimensional black holes is the existence of dynamical instabilities in 
the long wavelength, which is originally discovered by Gregory and Laflamme for the uni¬ 
form black string solution (UBS) whose horizon admits the isometry of S D ~ 3 x S 1 [1]. The 
Gregory-Laflamme (GL) instability happens if the circumference of the compact dimension 
L is in the same order of tq, which is the radius of S D ~ 3 . It is also known that the zero mode 
of the GL instability generates another static branch, called the non-uniform black string 
(NUBS), in which the translation symmetry in the compact dimension is broken. NUBS 
solutions are constructed perturbatively [2, 3] and then numerically up to D = 15 [4-8]. 

The phase of NUBS for the fixed L is known to admit a critical behavior with respect 
to the dimension [3] . In D < 13, NUBS becomes more massive than UBS at the GL critical 
point and has less entropy than UBS with the same mass. Therefore, the NUBS phase is 
unstable and cannot be the end point of the GL instability. On the other hand, for D > 13, 
NUBS becomes less massive and has more entropy, and then becomes stable. 

Recently, the authors and other people have shown that, in the large number limit of 
the dimension D —>• oo, the black hole dynamics can be solved analytically by the expanding 
in the inverse dimension 1/D [9-14]. These simplification can be explain by the fact that 
the large D limit of various black hole spacetimes universally admit the 2D string black 
hole structure, even in the existence of the cosmological constant and rotations [15]. 

The key feature of the large D limit is that the gravitational dynamics are divided into 
the two sectors each with the separated scales in terms of the black hole radius ro [12]: 

• Non-decoupling sector: the dynamics with the gradient large as ~ D/r o, which de¬ 
scribe the degrees of freedom of the gravitational wave propagating to the asymptotic 
region. 

• Decoupling sector: the dynamics with the smaller gradient (<C D/tq), confined within 
the short distance (~ tq/D) from the horizon, and decoupled from the asymptotic 
dynamics. 

The dynamics reflecting the characteristics of each black hole, including its instabilities, is 
described by the decoupling sector. As for the GL instability of the black string, the 1/D 
expansion revealed that the threshold wavenumber approaches /cql — y/D/r o at D —> oo, 
as suggested from the numerical result [3] . Its correction up to the fourth order in 1/D was 
also calculated and reproduced the numerical result, even for relatively lower dimensions 
D ~ 10 [9, 10, 14]. 

From the decoupling nature, one can expect that the dynamics of the small gradient 
degrees of freedom reduces to an effective theory living in the near-horizon geometry. The 
authors and collaborators have shown that the Einstein equation for the decoupling sector 
is separated to the radial direction with the large gradient d r ~ D/r o and the dynamics 
along the horizon with the smaller gradient (<C D/r q) [16]. Under this separation, the radial 
dependence can be easily integrated out by solving the ordinary differential equation, and 


- 2 - 


the remaining dynamics along the horizon is governed by a simple effective equation, which 
is expressed as the ‘soap-flim’ equation on the some constant surface of the radial coordinate 
placed at the short distance (~ tq/D) from the horizon. This formalism enables one to 
solve the Einstein equation even in the non-perturbative regime, provided the gradient 
along the horizon remains not so large (<C D/r o). The authors of [17] also studied the 
large D effective theory including the slow time dependence. 

In this paper, we study the effective equation for NUBS up to the next-to-next-to- 
leading order (NNLO) in the 1/D expansion, extending the leading order result in [16]. 
We also demonstrate the matching to the far-zone geometry and evaluate the asymptotic 
charges. With these results, the thermodynamic variables are calculated up to NNLO. 
The phase diagram including the NNLO correction admits the critical behavior around 
D* ~ 13.5. 1 Within the error of 0(D ~ 2 ), this value agrees with numerical observation [3, 8]. 

The construction of this paper is as follows. In section 2, we explain the setup and 
outline of the analysis. In section 3, the far-zone geometry is studied. In section 4, the near 
zone is solved by the 1/D expansion and the effective equation up to NNLO is derived. In 
section 5, we study the effective equation. In section 6, the thermodynamics of NUBS is 
studied. Section 7 gives the summary of this paper. The appendices include the details of 
some lengthy equations used in the main part. We attached a Mathematica file containing 
the detail of the NNLO near solution which are too lengthy to show here. The file also 
contains the series solution for the effective equation. 


2 Setup 

In this paper, we consider the D = n + 4 static and spherically symmetric spacetime with 
one compact dimension. Such spacetimes are described by the following metric 

ds 2 = -A(r, Z)dt 2 + B(r, Z)(dr 2 + dZ 2 ) + r 2 C(r, Z)dtt 2 n+1 , (2.1) 

where we assume A, B, C —i 1 at r —i oo for the asymptotic flatness. In this conformal 
ansatz, we solve the near-horizon geometry by the large D effective theory approach [16]. 
Because the gravity leaking from the near zone falls off by 1 /r n , we rather use 1/n as the 
small expansion parameter, instead of 1 /D itself. For the near-zone analysis, we introduce 
the another radial coordinate suitable to describe the near-horizon structure with the large 
radial gradient (~ n/ro), 


R = r n /rg. 


( 2 . 2 ) 


We also introduce the new compact coordinate to capture the variation in the scale of 

kch ~ V™/ r o, 



(2.3) 


1 In the literature, the critical dimension D* is usually defined as the highest dimension with the first 

order transition: D * = 13. In this paper, we adopt the definition of D* as the real number at which the 
order of the transition changes, as in the condensed matter physics. Because the large D analysis treats 
the number D just as a parameter, this is a rather natural definition. 
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This scaling is a natural choice if we want to keep the self-gravitating effect finite in the 
large n limit. 2 Assuming a bulge formed on the horizon with the wavelength A and radius 
deformed by <5ro ~ rg, the mass excess from the original horizon is roughly estimated by 


GSM ~ 


n+1 

2n 2 t-qA 

rp^) ' 


(2.4) 


Then, the extra gravitational potential made by this mass excess on its neighbor becomes 




GSM 

A n+1 





(2.5) 


which means that the gravitational interaction between bulges formed at intervals of A 
has the finite limit at n —>• oo only if A ~ rg/n 1 / 2 . For more gentle deformations like 
A ~ ro, the deformed horizon needs other supporting forces in the limit. For example, the 
authors and collaborators have demonstrated the negative cosmological constant allows 
several deformed solutions: droplets and cavities [16], and the rotation also allows bumpy 
Myers-Perry solutions [19]. As for the Kaluza-Klein spacetime, we have localized black 
holes other than black strings [7, 20], in which case the gravitational forces from distant 
black holes support the deformation. Though, such gravitational effect beyond the near 
zone will be rather contained in the non-decoupling sector with the large gradient. 

In our scaling (2.3), the static ansatz (2.1) allows the deformation only up to the mag¬ 
nitude of rg/n for the valid 1/n expansion [16]. However, this does not mean the analysis 
is perturbative. This magnitude of the deformation still produces the non-perturbative 
variation in the mass scale as (r o + Sro(z)/n) n ~ r n e Sr 0 (z)/r 0 _ \y e see ^jg non- 
perturbative degree of freedom follows a nonlinear effective equation. 


Matching Strategy For the evaluation of the asymptotic charges such as the mass and 
string tension, the author and collaborators proposed the general way of matching with 
the far zone through the quasi-local energy momentum tensor [16]. Instead, here we take 
the more primitive method. We will directly match the near and far solution in the overlap 
region 1 <C R <C e n . The strategy is as follows. The far solution is simply solved by the 
Newtonian approximation in the Minkowski background. The Newtonian potential <h(r, z) 
is expanded in the overlap region as 




( 2 . 6 ) 


where the function (j)(z) is identified with the mass deformation e 5r o(z)/uy Since 5ro(z) is 
only the functional degree of freedom in the near region, terms higher in InR/n will not 
provide any additional information. On the other hand, if we take r —>• oo, the Newtonian 
potential <h(r, z) admits the monopole behavior 


3>(r, z) 


arn 


(2.7) 


2 In the earlier work [10, 18], this scaling was also explained by the hydrodynamical picture. 
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The important fact is that the monopole a is obtained by averaging <f>(z ) over z, 

a = (</>) ~ ^ e ' 5ro(2) / r °y (2.8) 

Since the total mass (~ ar^L) will be the integration of the local mass density e Sro ^ z ^ ro , 
this is a natural consequence. Therefore, once the near deformation Sro(z) is solved, the 
asymptotic charges are determined. 

3 Far zone analysis 

In this section, we examine the far-zone geometry, to which the near solution in the 1 fn 
expansion will be connected. In the far zone r — ro ro/n, the gravitational potential 
~ rg/r” becomes much smaller than any power of n. Then, the far-zone geometry should 
be solved in the expansion of the Newtonian order r^/r”, instead of 1/n. 

We also study the boundary behavior of the far-zone geometry. The far zone includes 
two boundary regions: overlap and asymptotic region. In the overlap region, the far solution 
is further expanded by 1/n and matched with the near solution. In the asymptotic region, 
we merely take r —> oo without taking the large n limit to obtain the asymptotic charges. 

3.1 Newtonian approximation 

Let us perturb the metric (2.1) from the Minkowski background. 3 

ds 2 = -(1 + 5A)dt 2 + (1 + 8B)(dr 2 + dZ 2 ) + r 2 (l + 5C)dn 2 n+1 (3.1) 

where we write 

A = 1 + 5A, B = l + SB, C = 1 + 5C (3.2) 

and each correction is assumed to be in the Newtonian order 0(rg/r n ). Then, the equation 
reduces to the Laplace equation in the cylindrically symmetric space. 

^ d 2 + n + 1 d r + <9§^ <k(r, Z) = 0 (3.3) 

where d> is introduced by 

&4 = -$(r,Z). (3.4a) 

Other components reduce to 


AR r n 

SB = n + —T7*( r > Z )’ 

(n + 1 )r n n -f 1 

(3.4b) 

= ( 0 un + —T7*( r ’ Z ) + z )> 

(n- — l)r n n + 1 

(3.4c) 


3 The similar analysis is found in [21]. 
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where (3 is an integration constant, and \k(r, Z) follows 

(■dr + d 2 z )*(r, Z) = 0. 


(3.5) 


One can realize that ^(r, Z) is the degree of freedom of the conformal transformation. So 
we just set ^(r, Z) = 0. Assuming the periodicity Z —» Z + L, the general solution with 
the regular asymptotic behavior is given by the Fourier decomposition 

n oo n/2 

Z ) = Q ^l /2 K ™/ 2 (2njr/L) cos(2vr jZ/L) (3.6) 

3 = 1 

where K n / 2 (x) is the modified Bessel function of the second kind. The reflection symmetry 
Z —> — Z is also imposed. In the limit r —>• oo, only the first term can contribute to the 
asymptotic monopole, since K n / 2 (x) falls off exponentially at x —>• oo. 

Here we note that eq. (3.4) admits the consistency relations in the Newtonian order, 


5A + 2 5B + (n - 1 )5C = 0 


(3.7) 


and 

(n + 1 )5B + &A = —4/3/R. 

3.2 1/n expansion in the overlap region 

In the overlap region, we introduce the rescaled wavenumber 

27rro 


(3.8) 


k = 


y/nL' 


(3.9) 


With r = roR ly/?l , the Debye expansion formula leads to the expansion in terms of 1/n, 

2 u ‘ 2 (—2 j 2 k 2 + j 4 fe 4 ) In R + (j 4 /c 4 /2 — j 2 k 2 ) In 2 R 


K, 


l/2 (V^jkR 1/n ) = ^ (l - In R + 


n* 


(3.10) 


where fj(k ) is defined up to NNLO in 1/n as 


fj(k) — ( 1 + x— + 


x 1 + 


1 


6n 72 n 2 J y n \ejk 
j A k 4 — Aj 2 k 2 


n/2 


rk* 

e 2 


1 


4n 


n^ 


+ I -2 j 2 k 2 + 


5 j A k A 7fk 6 j 8 k 8 


12 


+ 


32 


(3.11) 


This expansion will be rather formal which breaks down for the large j. But we believe 
such high frequency components are negligible in the small gradient assumption. Then, 
the 1/n expansion of eq. (3.6) becomes 


<f> = — | a + ^ fj(k)aj{k) cos (jkz) + 0(ln R/n) 

3 =i 


(3.12) 
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where we also inserted Z = rpz/y/n. The leading behavior can be surned to a functional 
degree of freedom 

OO 

4>(z) = a + fj(k) a j{k ) cos (jkz). (3.13) 

3 = 1 

With this function, eq. (3.12) up to the relevant order can be rewritten in this form 

$ = I (fa) + ^\ z )— + $( 2 )(z)^ + 0(n~ 3 )] , (3.14) 

K \ n n z J 

where each is given by 

$ (1) (2) = 0"(z) + -(2 <j>"{z) + ^ 4 \z)), & 2 \z) = + \cf>^(z). (3.15) 

n Z 

We replaced the factor j 2 k 2 in eq. (3.10) with — d 2 . Therefore, the far solution (3.4) is 
determined only by < f>(z) and /3 in the overlap region. In the following, we see that /3 is 
related to the tension. 


3.3 Asymptotic charges 

Once 4>(z) and /3 are specified, the asymptotic charges are calculated as follows. In the 
asymptotic region r —>• oo, the leading behavior of 4> in eq. (3.6) is given by a, which is 
equal to the average of cf>(z) over z 


1 f L k i‘‘ 27T / k 

a = ((/)) = - J (j>(z)dZ = — J (p(z)dz 


The far solution (3.4) gives the monopole at r —>• oo, 


<5A = - 


arp §B = a - 4/3 rfi 


5C = 


8/3 


+ 


a 


77 + 1 r n 1 \n 2 — 1 n + 1 

Therefore, the ADM formula determines the mass and tension as [22] 


M = 


D n+ i rfiL 77(77 + 2) {(f>) + 4/3 


167rG 


77 + 1 


q- ttn+irp „ 

4v tG P ’ 


(3.16) 


(3.17) 


(3.18) 


where D n+ i = 2'K n / 2+l /T(n/2 + 1) is the volume of S n+1 . 


4 Near zone analysis 

In this section, we study the near-horizon geometry in the 1/n expansion, following the 
previous formulation [16] up to the one higher order. 
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4.1 Setup 

The near-horizon geometry can be separated into the radial sector with the large gradient 
(~ n/ro) and the other sector with the smaller gradient (<C n/ro), 


A r 2 

ds 2 = 1 V 2 (R, z) —77- + gu U dx fJ ’dx 1 '. 
n z 

Under this assumption, it is convenient to decompose the Einstein equation as 

K 2 - K^ u K p fi -R = 0 
S7 U K\ - = 0 

71 

j^d R K\ = -KK\ + R\ - -V^V.iV 
where R^ v is the intrinsic curvature for g /w and 

77 

K>\ = —g^d R g av . 

Our coordinate ansatz (2.1) leads to the intrinsic metric as 

r 2j, 2 

g IJV dx IJ dx 1 ' = —H(R, z)dt 2 + B( R, z) 0 “ —I- rQR 2 / n C(R, z)dti^ +ll 


n * 


with the conformal gauge 


N(R,z) = 


r 0 B(R,z) 

R 2 ^ 


(4.1) 


(4.2) 

(4.3) 

(4.4) 


(4.5) 


(4.6) 


(4.7) 


Since the far solution (3.4) has only the normalizable corrections (~ 1/R) in the overlap 
region, the boundary condition for R —>• 00 is simply 

A=1 + 0(R~ 1 ), 5 = 1 + 0(R~ 1 ), C = l + 0(R~ 1 ). (4.8) 

The large dimensionality of S'” -1 " 1 and the rescaling (2.3) also requires 

B = l + 0{n~ 1 ), C = l + 0{n~ 1 ) (4.9) 

for the valid 1/n expansion. 1 Hence, the metric components are expanded by 1/n as 

^*> = E 4 ^. b { r ,.-) = 1 + £c ( R,,) = 1 + E cm{Kz) 


k =0 


k =0 


n" 


k =0 


n 


k +1 


(4.10) 


4 The validity of the expansion itself requires only the constancy of the leading behavior in B and C [16]. 
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4.2 Leading order 

As was done in [16], the leading equation for the mean curvature K becomes 

—d R K( R, z) = -K( R, zf + (4.11) 

ro rg 

If we expand K as 


ro 4-^- n fc 
u k =o 


the leading solution iG 0 ' is given by 


(0) = R 2 + M{z) 2 
R 2 - M(z) 2 


(4.12) 


(4.13) 


where M(z) is an integration function denoting the horizon position. Solving eqs. (4.2) 
and (4.4), we obtain the leading solution 


A W = ( R M [ Z \ X, C® = 41n(l + 
VR + M(z)J V 


M(z) 

R 


and 


„foi - M(z)M"(z)) f M(z)\ 

M(z ) 2 { R J 


(4.14a) 


(4.14b) 


Adding to the condition (4.8), the regularity at R = M(z) is also imposed for B and C. 
The vector constraint (4.3) is trivially satisfied in this order. By setting M(z) = 1, this 
solution reproduces the uniform black string metric (A. 10). Here we note that M(z ) does 
not directly denote the deformation in the radius 5ro(z), but the variation of the mass scale 
M(z ) ~ e 5r o(z)/ra' 

Recalling the definition of R, the radius of the horizon embedding surface is given by 


rh(z) = r 0 M(z) l/n = r 0 exp 


— In M(z) . 
n 


(4.15) 


Since we assumed 5rh ~ 0(ra/n), the expansion is valid only if | \nM(z)\ <C n. 


Gauge choice As another gauge choice, one can set M(z) = 1 by the radial rescaling 
R —>• M(z)R and instead, use the integration function of C ^ as the degree of freedom [16], 

(4.16) 


C® = <?(,(*)+4 In ^1 + 


This choice makes the radial coordinate fit the equipotential surface, which leads to the 
simpler near horizon analysis. 5 In this ‘equipotential’ gauge, the near deformation function 
enters in the non-normalizable behavior in the overlap region, C'( 0) = Cq(z) + 0(R _1 ). So, 

5 In the current system, the similar coordinate was used by Harmark and Obers [23]. 
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in exchange for the simplified near-zone analysis, the matching with the far-zone geometry 
requires a nontrivial coordinate transformation, involving both R and z. 

Instead, the current embedding condition (4.8) makes the near coordinate rather fit 
the asymptotic geometry. The deformation in the non-normalizable behavior Cq{z) is now 
absorbed to the horizon shift M(z). This choice is convenient to keep contact with the 
asymptotic coordinate, and in which the matching between the near and far zone will not 
require any further coordinate transformation other than R = r n /rQ. 6 


4.3 Next-to-leading order 

The solution at NLO is similarly calculated. The horizon position is fixed to R = M ( z ) by 
the renormalization of M(z). In the appendix. B, we show the detail of A^\ B ^ and C^\ 
At this order, the vector constraint (4.3) admits the first nontrivial condition for M(z), 


0 = M®\z) + M'(z) 


2 M'(z)M"(z) M'(z) 3 

M{z) + M(z) 2 ' 


(4.17) 


This is equivalent to the differentiated version of the effective equation obtained in [16] . 7 As 
pointed out in [16], eq (4.17) is the system equivalent to the undamped Toda oscillator [24]. 


Effective equation from the matching Due to the assumption dr^(z) ~ ro/n, the 
vector constraint (4.3) provides only the one lower order effective equation. This situation 
can be remedied by the use of the far-zone information. We can also derive the leading order 
effective equation only from the leading solution (4.14), just by using the far consistency 
relation (3.8). The leading solution (4.14) expanded up to 0(R _1 ) is matched with the far 
Newtonian correction (3.4) as 


<5 A = 


4 M(z) 4(M'(z) 2 - M(z)M"(z)) 4 M(z) 

R ’ “ nM{z )R ’ C “ nR 


Substituting this into eq. (3.8), we obtain the leading order effective equation 


M(z) 


M'{z) 2 
M(z) 


+ M"(z) = 0. 


(4.18) 


(4.19) 


This is just the integrated form of eq. (4.17). 

We note that this derivation is due to the fact that the far-zone equation is solved 
explicitly, which depends on the asymptotic structure. More general asymptotics will not 
allow such simplification. 


4.4 Next-to-next-to-leading order 

Finally, we obtained the solution up to the next-to-next-to-leading order (NNLO). Because 
of its lengthy form, the detail of the NNLO solution is given in the attached Mathematica 

6 This only works if the embedding has the trivial leading behavior: r = ro + Ofn -1 ). For the nontrivial 
embedding like r = IZ(z) + Dfn" 1 ) as considered in [16], the use of M(z) will not simplify the matching. 

7 One can find that M(z) corresponds to e 2V ^ in [16]. 
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file. To this order, the vector constraint (4.3) gives 1 


. .. . 2 M'(z)M"(z) M'(z) 3 

1 j ^ 1 1 M{z) M{z ) 2 


1 

n 


(5 - 2C(2) + 4 In 2) 


, (2C(2) -2-81n2 )M\z) 3 

+ M(ip + 


M"(z) 2 M'(z) M'{zf 
M(z) 2 + M(z) 4 

(81n2-2C(2))M 


Af(«) 


^ + (1 + 41n2 + 21nM(z))M , (z) 

'( 2 ) (4C(2)-10-81n2)M'(z) 3 

+ Wzj 3 


M"(z) 

(4.20) 


where the derivatives higher than M"(z) in 0(n -1 ) are eliminated by using this equation 
iteratively. This effective equation determines the horizon deformation only up to NLO. 
We will see that the one higher order equation can be obtained from the matching. 

4.5 Matching 

Now, we identify the Newtonian correction in the far zone (3.2) with the near solution up 
to NNLO. Expanding the near solution up to 0(R _1 ), we obtain 

2 a(£) i 

y2 — = 1 - p {-Ao(z) + Ai{z) In R/n + A 2 {z)(\n R/n ) 2 ) + 0{ R~ 2 ), (4.21a) 
i=0 

X + = l + ^E{Bo(z)+B 1 (z)lnR/n + B 2 (z)(lnR/n) 2 ) + e>(R“ 2 ),(4.21b) 

n z ' rr nK 

i=0 

1 + - Y" ——j— = 1 + {Cq{z) + Ci(z) In R/n + C 2 (z)(In R/n) 2 ) + 0(R~ 2 ), (4.21c) 
n n l n K x 7 

»=o 

where each coefficient consists of M(z ) and its derivatives as in eq. (4.18). We show only 
the detail of Ao in the appendix C. The Newtonian correction 5A, 5B, 5C in the far 
zone (3.4) is matched with 0(R _1 ) terms in the above equation. As in the leading order, 
eq. (3.8) for 0((lnR)°/R) gives the effective equation up to NNLO 

Ao(z) - (n + 1 )B 0 (z)/n = 4/3. (4.22) 

The detail of this equation is also shown in the appendix C. The vector constraint (4.20) is 
equivalent to the derivative of this equation up to 0(n -1 ). The match of 5A also determines 
cf)(z) as 


<j>(z) = Ao(z). (4.23) 

Therefore, we have completely determined the far-zone geometry up to NNLO in 1/n. 

8 The relation between this equation and the soap-film equation in [16] is not trivial, because our radial 
coordinate is not normal to the horizon. We need the transformation to the normal coordinate to find the 
relation. 
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Consistency check Eq. (4.21) also contains the power of In R/n. Since we have no more 
undetermined degree of freedom, the constraints (3.7) and (3.8) should be trivially satisfied 
for such terms. Actually, we found that our near solution satisfies up to the relevant order 

A k - 2 B k /n - (n - 1 )C k /n = 0 [k = 0,1, 2) (4.24) 

and 

A\ — (n + l)B\/n = 0, A 2 — (n + 1)32/n = 0. (4-25) 

Using eqs. (4.20) and (4.23), the remaining consistency (3.15) is also satisfied. 


5 Analysis of the effective equation 

In this section, we study the non-uniform black string solution by solving the NNLO 
effective equation (4.22). 


5.1 Integrated form 

After an integration, eq. (4.22) is arranged to the following form 


0 = 

+ 


M'(z) 2 , , . In 2 M(z) In 3 M(z) 

W — a 0 — a! In M(z) — a 2 -—^ - a 3 - 


2 M(z) 2 

P 


n 


M(z ) 


bo,o + bo.i- 


ln M (z) 


n 


+ bi 


In 2 M(z) 
0,2 - 5 — 


+ 


P 


nM(z) 


bi,o + bii - 


In M(z)\ b 2 ,o ft 


n 


+ 


n 2 M(z) 3 

(5.1) 


where /3 is the rescaled value of tension in eq. (3.18). Each coefficient is given by 


3o — 1 + CL -f- 


1 + 4a + 2a 2 + 4 In 2 1 


2 n 




4 —y ( a + 2a 2 + -a 3 + 2(2a + 1) In 2 ) , 


1 4(a + l) 2(1 + 4a + 2a 2 + 4 In 2) „ 8(a + l) 14 

a i — — 1 - 5 - , a 2 — 3 H-, 33 — — 


n 


n z 


n 


(5.2a) 


3 ’ 


bo.o = 1 + 


3 + 4 In 2 - C(2) - 2a(C(2) - 2) 


n 


+ 


10 n 2 


(30 + 80 In 2 2 - 80(C(2) - 2) In 2 - 70C(2) + 49C(2) 2 + 4a 2 (20 - 25C(2) + 7C(2) 2 ) 


- 15C(3) - 2a (-60 + 40(C(2) - 2) In 2 + 80C(2) - 42C(2) 2 + 15C(3))) , 
b 0 ,i = 2C(2) - 6 


(5.2b) 


1 


n 


42 


4a, 


+ - -18 + 8(C(2) - 3) In 2 + 18C(2) - — C(2) 2 + 3C(3) - —(30 - 30C(2) + 7C(2) 2 ) , 


14 


bo. 2 = 18 - 14C(2) + ^C(2) 2 , 
5 
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5 1/ 46 

b, 0 =-2 In 2 + C(2) - - 10 + 8 In 2 2 + (26 - 20C(2)) In 2 - 13C(2) + — C(2) 2 

2 n \ 5 

+ a ^15 - 8 In 2 2 - 16C(2) + y C(2) 2 - 4(4C(2) - 3) In 2 + 2C(3)^ ^ , ( 5 - 2c ) 

bi i = 20 - 8 In 2 2 - 16(C(2) - 1) In 2 - 18C(2) + — C(2) 2 + 2£(3) 

5 

and 

24 3 

b 2 , 0 = 7 - 12(C(2) - 1) In 2 - 6C(2) + —C(2) 2 + -C(3). (5.2d) 

5 2 

where a is an integration constant, whose definition is tuned for the later convenience. This 
is one of our main results. Here In 2 in the coefficients is just originating from the current 
conformal coordinate, and does not affect the physical quantities. 

Eq. (5.1) is invariant under the rescaling, 

M(z) Mi(z\) = e^M(e~^ n zi), /5 —>• /?i = e^/3, a —>• a\ = a + £. (5-3) 

The definition of a is arranged so that this transformation works linearly on a. Then, the 
constant a is found to be the scaling degree of freedom of the spacetime, which can be set 
to an arbitrary value. The scaling at a = 0 is also fixed so that (3 = 1 just admits the 
uniform solution. 


5.2 Leading order equation 


Let us start with the leading order equation 


In M + 


M' 2 
2 M 2 




(5.4) 


For the time being, we put aside the scaling a by setting a = 0. If we define <p(z) = y/M[z ), 
then eq. (5.4) reduces to the potential problem of the classical dynamics, 

\p’ 2 = V{(p) = J(v? 2 hV - <^ 2 ) (5.5) 


where the minus of the rescaled tension — /? plays the roll of the energy. The minimum 
energy state with ft = 1 is given by the uniform solution <p(z) = 1. For the positive 
tension 0 < f3 < 1, this potential also admits the oscillatory solution, which corresponds 
to the non-uniform black string (NUBS). Since the expansion breaks down at M(z) = 0, 
the solution for the negative tension /3 < 0 is not the physical branch. The solution also 
collapse to M(z) = 0 for the zero tension f3 = 0. However, the zero tension solution can be 
identified with localized black holes, which will be discussed in the later section. 

The maximum and minimum value for NUBS are given by 


^7{ maXjm i n j (/?) 


^{ 0 ,- 1 } (-/5/e) 


(5.6) 


where Wm-i}^) denote the upper and lower branches of the Lambert W function, defined 
as the inverse function for z(w) = we w . In the limit /3 —>• +0, we have 


M max ~ e, M min ~/3/ln(l//3). 


(5.7) 


- 13 - 





Since the expansion is only valid for | \nM(z)\ -C n, the range of the parameters is restricted 
to /3 e~ n . 


Non-uniformity parameter An useful dimensionless measure for the deformation is 
the non-uniformity parameter [2] 


A =2 


1 /fl„ 


R n 


- 1 


(5.8) 


where f? max and -R m in is the maximum and minimum radius of S n+1 . This is expressed by 
the maximum and minimum value of M(z ), 


A = 


1 ( M max (f3) 1 / r ‘ 


- 1 


(5.9) 


M min (^)Vn 

Using eq. (5.6), the non-uniformity parameter is estimated for small f3 and large n as 

(5.10) 




where we assumed f3 e n . Thus, the l/D expansion covers the deformation only with 
the small non-uniformity A <C 0(1). 


5.3 Higher order corrections 

Now, we study the effective equation including higher order corrections. For a = 0, eq. (5.1) 
is arranged to gives the UBS solution for ft = 1 as in the leading analysis, 

Mubs = 1 - + ^2 (C(2) 2 - 4C(3)) • (5.11) 

The correction reproduces eq. (A. 11), if we identify = ro(-MuBs) 1 ■ In the following 
analysis, we introduce the scale invariant variables M(z) and b dehned by 

M{z) = e a M(e~ a/n z), (3 = e a b. (5.12) 

For a = 0, M(z) and b coincide with M(z) and /3. The boundary values for M( 2 :) is given 
by 

+ I2S?( 5(11_721112 + 48 ‘“ 22) 

+ (65 + 120(—5 + 4C(2)) In 2 + 720 In 2 2 - 240C(2) - 84C(2) 2 + 180C(3))W {0 ,_i} (- b/e) 

+ 2(5 + 120(—1 + 2C(2)) In 2 + 240 In 2 2 - 120C(2) - 12((2) 2 - 30C(3))fU {0 _i } (-6/e) 2 ) . 

(5.13) 


1 + i- (4In 2 - 1 + (2C(2) - 1 + 4In 2) W" {0> _ 1} (6/e)) 


One can see that, by the scaling property, the non-uniformity parameter (5.9) depends only 
on the scale invariant parameter b, 


x= l( M max (fr) 1 / n _ \ 

2 ^ M min (6)Vn ) ■ 


(5.14) 
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Near uniform solution Since we cannot take 6 to be as small as e~ n , it suffices to 
consider the expansion around UBS, 6=1. Let us go back to the leading equation (5.4). 
Expanding with 6=1 — e 2 /2 and M(z) = 1 + <5M(z), eq. (5.4) admits the small oscillation 
in 0(e) 

6 M(z) = ecos(z) (5.15) 

where we assumed the reflection symmetry at z = 0, M(z) = M(— z). Since the proper 
length is given by Z = r^z/y/n, this reproduces the leading behavior of the GL threshold 
wavenumber /cgl — y/n/rQ. Due to the nonlinearity in eq. (5.4), the period is also corrected 
for the absence of the secular terms as zsin(z), 

cos( 2 :) —>• cos (k(e)z). (5.16) 

The first correction to k(e) is obtained at 0(e 2 ) 

k(e) = l-^e 2 + 0(e 4 ), (5.17) 

which is determined by the solution up to 0(e 3 ), 

e 2 e 3 / 7 \ 

M(z) = 1 + ecos (k(e)z) H-cos(2 k(e)z) - ( - cos (k(e)z) — cos(3 k(e)z) ] + 0(e 4 ). 

6 96 \ 3 J 

(5.18) 


The inclusion of higher order corrections in eq. (5.1) is rather straightforward. We just 
expand M by 1/n as M = Mo + Mi/n + M 2 /n 2 , and the avoidance of the secular terms also 
adds the higher order corrections to the wavelength k(e) = ko(e) + k\(e)/n + k 2 (e)/n 2 . If 
we continue the analysis, the resulting solution can be written in the Fourier series 

OO 

M(z) = cos(jk(e)z). (5.19) 

l=o 


where the j-th harmonics begins from 0(e J ). We show only the first few terms in the 
appendix. D. Further terms up to 0(e 7 ) are given in the Mathematica file attached to this 
paper. 

The wavenumber k(e) determined by the expansion up to 0(e 7 ) is given by 


1 


k(e) = 4-' ln ( 1 - — + 
w 1 2n 8n 2 / 


1 - 11 -^+ 7 
n n 

379 




9 


+ 


6971 

2488320 


2 7 24 V 768 128n 256n2 

23 \ 


-e b + 0(e 8 ) 

18432n 41472n 2 ) [ ’ 


(5.20) 


In the UBS limit e —>• 0, k(e ) reproduces the higher order corrections to /cgl shown by 
the linear analysis [9, 10, 14], except the factor 4~ 1 / n . Actually, we will see that the mass 
length of the corresponding UBS is given by po = 4}^ n e a ^ n rQ. Then, if we recover the 
scaling (5.12), the physical wavenumber at e = 0 is given by 


y/nk{ 0) 

e a / n r 0 


kch 


y/n 

Po 




(5.21) 
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In (M) In (/W) 




Figure 1. The convergence of In M(z) for n = 10, e = 1.2 with respect to e (left) and 1/n (right). 
Only the half period is plotted. In the left, each curve corresponds to the NNLO solution up to 
0(e 3 ), 0(e 5 ) and 0(e 7 ), respectively. In the right, each curve is the series solution up to 0(e 7 ) 
including the LO, NLO and NNLO correction. While the series in e converges, the expansion in 
1/n is not so good around the neck In M ~ —4. 


The series solution in eq. (5.19) itself is convergent even around e ~ \/2 where the 
rescaled tension b = 1 — e 2 /2 reaches zero, regardless of the convergence in the 1/n ex¬ 
pansion (figure. 1). So, in the following analysis, we use this series solution, instead of 
solving eq. (5.1) numerically. Since the parameter a just involves the scaling, the shape 
of the deformation is governed only by e. Figure 2 shows that as e approaches to \/2, 
the deformation tends to have the broader bulge and narrower neck as seen by the nu¬ 
merical calculations. Analytically, this can be shown by estimating the inflection point 
<9 2 In M(zo) = 0, 

k{ e ) z o - | = e Q + 7^(1 + 2C(2) - 4In2) 

+ 2 fj 2 + 3 C( 2 ) + 5 C(2) 2 + 10(1 - 2C(2)) In 2 - 12 In 2 2 + 7C(3)) ) + 0(e 3 ) 
^(l + - + —y 3+ 0(e 3 ). (5.22) 


6 Thermodynamics 

In this section, we study the thermodynamic properties of NUBS in the 1/n expansion. 

6.1 Variables 

First, we compute the thermodynamic variables from the solution (5.19). The mass and 
tension are derived from eq. (3.18). Using eqs. (4.23) and (C.l), the average of 4>(z) can 
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In (M) 



Figure 2. The deformation for n = 10 in a half period. The curves are for e = 0.4,0.6,0.8,1 
and 1.2, each with the non-uniformity A = 0.0566,0.0912,0.135,0.197 and 0.316, respectively. The 
corresponding non-uniformities are calculated by eq. (5.14) with b = 1 — e 2 /2. Each dot on curves 
represents the point where the sign of d 2 In M( 2 ;) changes. 


be expressed by the average of M(^) 


<«/>) = (A) = 4e° 


<M> + 


(M) + 6(C(2)-1) 


n 


+T2 C(2) - xC(2) 2 + 2C(3) - 1 - 2C(2) + —C(2) 




7 


10 



( 6 . 1 ) 


where we eliminated lnM(z) and higher derivatives of M(z) by using eq. (5.1) and inte¬ 
grating by part. Using the series solution (5.19), the above averages can be computed 
as 


(M) = /i 0 (e), 



( 6 . 2 ) 


For the latter, we also used /xi(e), yr* 2 (e), A* 3 (e) and eq. (5.20). The form of /xo(e), /ii(e), // 2 (e) 
and // 3 (e) is given in the appendix. D. In the UBS limit e —>• 0, we obtain the mass of UBS 
on the GL threshold point 


A4gl = lim M 

e —>0 


(n + l)U n+ iLe a ro 
4 t tG 


(6.3) 


Comparing with eq. (A.2), the length scale of the UBS mass on the threshold becomes 


Po = (4 e“) 1 /n r 0 . 


(6.4) 


On the other hand, since the horizon position rh(z) = roM(z) 1//ri varies in z, the horizon 
area is given by 


A = r£ +1 fl 


n +1 


r b H {z) i / 2 c j 

Jo 



n±i M(z) 2 / n (M'(z)y 


nM(z) 


M(z) 1+1/n dz (6.5) 
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where Bh(z) and Ch{z) are the metric components on the horizon. Using eq. (5.1) and 
integrating by part, the area just reduces to 


A = A 1 / n e^n n+1 r^ +1 L(^). ( 6 . 6 ) 

where we used the form of eq. (6.1). The surface gravity is also computed by using eq. (5.1) 

(6.7) 


n 


k = 


2 (4e a ) 1 / n ro 

The Smarr’s formula just follows from eq. (3.18) and the above expressions for A and n, 


(n + 1 )M - TL _ nfln+irgL _ kA 
n + 2 167 tG 87 tG 

In fact, the scaling factor can be factored out from these variables 

JVn(4 1/n ro) n+1 «±i 0 (n + 2)n^ 0 (e) + 1 - e 2 /2 


( 6 . 8 ) 


M = 
T = 


I 67 tG 
^n+irft 
4t tG 


n + 1 


Lo(e), 


e a (l - y) , A = O n+ 1 (4 1 /- ro ) ri+ 2 e^>o(e)io(6), 


where we introduced the scale invariant quantities 

■ / N _ (4>) T ( \ - L 27F 

^ J _ 4e“’ ° U _ (4e a y/ n r 0 4 
Substituting eq. (6.2) into eq. (6.1), ^>o(e) up to 0(e 6 ) is given by 

e 4 lie 6 1 /e 2 e 4 79e 6 \ 1 / e 4 29e 6 \ 


4(e) = 1 - — - 


+ 


+ — + 


96 3456 ' n V 2 ' 12 ' 3456 ) 

First law Now, we check the first law for NUBS 

cLM = ~^dA + TdL. 

87rG 


n 2 V 48 + 3456 ) 


(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


( 6 . 12 ) 


The first law for the scaling a is easily confirmed. For the variation with e, the first law 
reduces to the condition 


- e - + (Mt) - ! + = 

de n de 


(6.13) 


From eqs. (6.11) and (5.20), this condition can be confirmed at least for up to 0(e 6 ) and 
0(n~ 2 ). 


6.2 Circumference 


So far, we have considered the scaling and deformation separately. However, we have to 
vary the scaling a together with e to obtain the physical phase diagram. If we recover the 
scaling, the circumference of the extra dimension is given by 


27re a / Tl ro 

y/nk{e) 


(6.14) 
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Therefore, the circumference also changes as the solution gets deformed. To obtain the 
sequence of NUBS with a fixed L, the scaling a should be tuned so that the running of e 
is canceled. Since the running in k(e) becomes comparable to e a ( n only if e 2 ~ n _1 , we 
introduce a rescaled non-uniformity parameter 

e 2 = -, (6.15) 

n 

where we assume s ~ 0(1). With this assumption, the variation in eq. (5.20) is cancelled 
by setting 


a ““l2 + n 


/3s 

23s 2 A 

1 | 

(7s 

11s 2 997s 3 

V 4 

576 ) 

n 2 ' 

V 12 

32 + 38880 


(6.16) 


where the scaling at s = 0 is absorbed in tq. Therefore, the physical phase of NUBS is 
expressed in terms of an 0(1) non-uniformity parameter s. 

6.3 Phase diagram 

Now, we investigate the phase diagram for the constant L. With the assumption e 2 = 2 s/n 
and the scaling in eq. (6.16), L is now fixed to 

2vr 


L = Lql — 


kch 


(6.17) 


where &gl is given by eq. (5.21). From eqs. (6.11) and (6.16), eqs. (6.9) and (6.7) gives 
M _. /19 (_ (_ _5_\ 3s /_23_ _7_\ 17291s 3 , 529s 4 \ 


M G l 

T 

7gl 


= e~ s / 12 


1+1 + 


+ 


9n J An V 576 12n J n 311040n 2 663552 n 2 ) ' 

(6.18a) 


= e ~ s / 12 ( 1 - ( 1 + 


7 


f 23 


+ 


1 


3n J 4n \576 8 nJn 




4871s 3 529s 

311040n 2 + 663552n 2 ) 




A = —s /12 ( x ( 2 7 '\ 1 £ _ 127s 3 529s 4 

1 3 6 n) n V 576 64 n) n 2430n 2 663552n 2 ) 

9041s 3 


Agl 
k 

kgl 

where the value at the GL point is given by 

_ (n + l)fi n +iPo^GL 


A 9 

7 \ 

s ( 

1-h 



\ n 

n 2 ) 

1 12n \ 


A4 G l 


1 1 ):: (i 

„n+1 


+ 


Tgl 


n 


311040n 3 ’ 

M gl 

(n + 1)L G l’ 


(6.18b) 

(6.18c) 

(6.18d) 


(6.19) 


AqL — ^n+lPo ^GL, «GL = -—■ 

4/50 


In figure 3, we compare our result with and the numerical result in [8] at several 
dimensions. The corresponding non-uniformity parameter A is determined by eq. (5.14) 
with b = 1 — s/n. The variables are plotted in the unit of Tgl 5 


M 


(^GL) n+1 

A 

(L G L ) n+2 


(n + l)Q n+1 (kGL] n f M ^ 


= n 


n +1 


fcGlA 

2vr J 


2vr 

n+1 


T 


Hn+l (^GL 


A 


\MglJ ’ (Lgl ) 1 

TITT 


167T V 27T 


T 


Tgl), 


4 


Agl) ’ 


kL G l = 


k G L V K GL 


( 6 . 20 ) 
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Figure 3. Phase diagram of NUBS, and comparison with the numerical construction of [8] (the 
data for the tension, not included in [8], has been provided by Pau Figueras). The normalized 
surface gravity k = n/n (top left), horizon area A (top right), mass A4 (bottom left) and tension 
T (bottom right) are plotted against the non-uniformity parameter A. Each quantity is normalized 
in the unit of the extra dimension size L. The dotted, dashed, dot-dashed and solid curves are the 
variables for n = 7, 8,9,10 (D = 11,12,13,14), respectively. 


where we write A;gl = Po^GL and set G = 1. The formula in eq. (6.18) reproduces the 
numerical result within the error of 0(n ~ 3 ) for the small non-uniformity. For the large 
non-uniformity, the 1/n expansion breaks down and one can no longer expect the good 
accuracy. 

In figure 3, the tension seems to have the relatively large error compared to other vari¬ 
ables. Actually, the Smarr formula (6.8) implies that the tension has the error of one lower 
order than other variables. This is because the tension is sensitive to the inhomogeneity of 
the horizon, especially the neck part of the horizon which has relatively bad convergence 
(figure 1). On the other hand, the mass and area, which are the averaged quantities, will 
not be so sensitive to the neck geometry. The surface gravity, which is defined uniformly 
on the horizon, also will not have such sensitivity. 


6.3.1 Localized black holes 

We should note that our formalism also includes localized black holes as the zero tension 
solution. If the solution has the point M(z ) = 0, the regularity in eq. (4.14b) requires 



M(z) ) M ( z ) =0 


( 6 . 21 ) 
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From the effective equation (4.19), this is the case of the zero tension (3 = 0. For (3 = 0, 
eq. (5.1) admits an analytic solution 


M(z) = e ~ z2/2 



z 6 (—16 + 3z 2 ) 

96+ 


( 6 . 22 ) 


where we omit the scaling for simplicity. This solution is not periodic nor convergent at the 
large 2 , so one can think this unphysical. But, this can be interpreted as the near-equatorial 
part of the spherical solution. 

Let us consider D = n + 4 Schwarzshild spacetimes in the isotropic coordinate 


ds~ = — 


n +1 


Po +1 


P n+1 + Po 


n +1 


dt + I 1 + 


n+l' 


4 

n+l 


(dp 2 + p 2 dd 2 + p 2 sin 2 9dQ 2 


n+l 1 


(6.23) 


We consider the transformation to the cylindrical coordinate 

r = p sin#, roz = \/npcos 9, (6.24) 


where 2 coordinate describes the region around the equatorial plane 9 ~ n/2. Then, we 
obtain 


ds 2 ~ — 


R-e-* 2 / 2 ' 
R + e-^/2 


dt + I 1 + 


0 -z 2 /2 \ n+l 


R 


2 , 7 0 


rZdz 2 


dr z + 


n 


+ r 2 1 + 


= — 2 2 /2 \ n+l 


R 


d^l+i 


(6.25) 


where we set po = ro and the following formula is used 


sin n+1 9 = [ 1 + 


9 \ 

■yZ \ 2 


n 


R 2/ r 


~ e 2 + 0(n ) 


(6.26) 


This is identical to what we have obtained in eqs. (4.14a) and (4.14b) with the leading 
order of eq. (6.22). 

As for the higher order, though it is difficult to find the coordinate transformation, 
Eq. (6.22) reproduces the mass spectrum of localized black holes, correctly up to NNLO. 9 
Since we have (3 = 0, substituting eq. (6.1) into eq. (3.18) gives 


M = ( w + 2 ) n "+iPo L (m) 

47 tG 

= (n + 2)1W£ +1 M{z) dz_ 

4v tG J~ y/n 

_ (n + 2 )tt n +iPo +1 [2 t[ ( _ 3_ 25 \ 

47 tG V n \ 4n 32n 2 / 


9 We thank Roberto Emparan for suggesting this possibility. 


(6.27) 
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where we identified L(M) with the integration over Z = r^z/y/n in the middle line. Using 
the following identity 


3 

V n \ 4n 32 n 2 J ' 


(6.28) 


eq. (6.27) reproduces the mass of D = n + 4 Schwarzschild black holes in the isotropic 
coordinate (6.23) up to NNLO, despite the solution (6.22) itself breaks down at the large 
z. As discussed in [16], this is because the contribution from the equatorial region becomes 
dominant for the energy density in the large D limit. 

However, our construction fails to capture the thermodynamics and deformation in the 
localized black hole phase. Eq. (6.27) is only the mass of the spherical black hole without 
the deformation. If we compere the mass with the NUBS and UBS phase for the same 
L, any size of the localized black hole has the zero mass, because of the infinite period. 
Then, in the large D limit, the localized black hole phase is degenerate to a single point 
of the phase diagram: (A4 = 0,T = 0). This is because the deformation effect is non- 
perturbative in 1/n for localized black holes, at least, in the small gradient assumption. 
The deformation of the localized black hole with the radius po and extra dimension size L 
is supported by the gravitational potential from itself at a distance of a half period L/ 2. 
This gravitational effect has the magnitude of (2po/L) D , which is negligible in the 1/D 
expansion, unless the black hole is very large L/2 — po Po/D- 


6.4 Critical dimension 


The phase diagram of NUBS admits a critical behavior between D = 13 and D = 14, across 
which the order of the transition changes [3] . This means that NUBS has the greater mass 
than UBS on the GL point for D < 13, and the lower mass for D > 13. In eq. (6.18a), the 
gradient near the bifurcation point becomes 


A4ql V n n 2 J 12’ 

in which the gradient changes the sign at 


(6.29) 


n 2 — 9n — 5 = 0. 


(6.30) 


This implies the existence of the critical dimension at n* ~ 9.5 (77* ~ 13.5), which agrees 
with the Sorkin’s numerical result [3], within the error of 0(n~ 2 ). Here we note that the 
n = 9 case admits a marginal behavior at NLO, in which case the NNLO correction plays 
a crucial role. 10 In figure 4, we show the critical behavior in the phase diagram, plotted in 
terms of the normalized mass and relative tension. The relative tension is a dimensionless 
quantity defined by 

(n+l)LT _ (n+l) 2 (l-e 2 /2) 

' M (n + 2)n(t>o(e) + 1 — e 2 /2 ’ 

10 A subtlety is that the leading and sub-leading terms become comparable around the critical dimension. 
One may suspect the validity of the expansion. But we believe this only happens at the leading order, since 
the subsequent correction seems convergent. 


- 22 - 








which is normalized so that the UBS limit gives r = 1. Substituting e 2 = 2 s/n, we obtain 



1 

1 H— 
n 




7 s 3 

432ra 3 ’ 


(6.32) 


where the corrections can be determined up to 0(n~ 3 ), because 0{n ~ 3 ) correction in 
eq. (6.31) should start with e 2 for the normalization r(e = 0) = 1. Substituting this 
expression into eq. (6.29), we obtain the gradient in the phase diagram, 


M 

-A/Igl 



, 10 6 

1-1-2 

n n z 



T )- 


(6.33) 


The condition n 2 — lOn + 6 = 0 also gives n* ~ 9.4. Since the gradient in eq. (6.32) does 
not affect the sign, this value is identical to n* ~ 9.5 up to the error in 0(n~ 2 ). 


04 



Figure 4. Critical behavior of NUBS near the bifurcation point, in which the normalized mass is 
plotted against the relative tension. The thick, dot-dashed and dashed curves are for n = 8,9,10 
(D = 12,13,14), respectively. We plotted the results up to the NLO (gray) and NNLO (black) 
corrections, in comparison with the numerical data for n = 8 (circle), n = 9 (square) and n = 10 
(triangle) in [8] (the data for the relative tension, not included in [8], is produced with the tension 
data provided by Pau Figueras). At NNLO, the mass seems to grow by the effect of the self¬ 
interaction, which makes the n = 9 phase bellow the critical dimension. 


The horizon area is also known to admit the critical behavior at the same dimension 
as the mass. Above the critical dimension, the horizon area of NUBS becomes greater than 
that of UBS with the same mass, which means that NUBS is more favored than UBS in 
the microcanonical ensemble. Below the critical dimension, NUBS becomes the unstable 
phase. However, we cannot observe the difference in the area between NUBS and UBS of 
the same mass to this order, 


■A/Axibs = 1 + 0(n 3 ) 


(6.34) 












where the horizon area of the UBS with the same mass is given by 


Ajbs 


^n+l-^GLPo +1 


n +1 , 


(n + 2)n</> 0 (e) + 1 
(n + l) 2 



(6.35) 


This is expected since the first law guarantees the cancelation of the linear order depen¬ 
dence [2, 3]. In this case, both A and A.ubs are the function of e 2 , and then the second 
order becomes e 4 . In addition, the coefficient will be at most since there is no en¬ 

tropy difference in the large D limit [10]. Then, the leading behavior of the area difference 
should appear at much higher order than the mass 


A/Abbs — 1 ~ 0(e 4 /n) ~ 0(n 3 ). 


(6.36) 


The numerical calculation actually shows -4 ./.Tubs — 1 — 10~ 3 ~ 1CU 4 for the small non¬ 
uniformity in D = 11 to 14 (See figure 7 in [8]). 

As for the surface gravity in eq. (6.18d), though its thermodynamical interpretation is 
unclear, we see that the gradient at the GL point changes at 


n 2 — 9n + 7 = 0, 


(6.37) 


which gives a smaller critical value n* ~ 8.1 (D* ~ 12.1). Then, the NUBS phase becomes 
hotter in D > 12 and cooler in D < 12. This effect is also observed by the numerical 
calculation [3, 8]. 


7 Summary 

In the large D limit, the near-horizon dynamics with the small gradient (<C D/tq) along 
the horizon reduces to the nonlinear effective theory for some collective degrees of freedom 
of the horizon [16, 17]. This effective approach will be a powerful tool for the various 
nonlinear problems in the higher dimensional gravity. 

In this paper, we applied the large D effective theory approach to the study of the 
non-uniform black strings (NUBS) bifurcating from the uniform black strings (UBS) on 
the threshold of the Gregory-Laflamme (GL) instability. Extending the leading order 
analysis [16], we solved the near-horizon geometry up to the next-to-next-to-leading order 
(NNLO) in the 1/D expansion. We also analyzed the far-zone geometry by the Newtonian 
approximation and derived the matching condition between the two zones. As the result, 
we obtained the effective equation of the deformed horizon up to NNLO. We also obtained 
the thermodynamic variables up to NNLO from the matching result. The phase diagram 
produced by the large D analysis agrees with the numerical results for the small non¬ 
uniformity within the expected precision. We also show that our construction contains 
localized black holes which have the degenerate phase in the 1/D expansion. The phase 
diagram admits the critical behavior in the mass around D * ~ 13.5 as shown numerically [3, 
8]. We cannot find the area difference between NUBS and UBS of the same mass in 
this order. It will be straightforward (though much harder) to obtain much higher order 
corrections, in which the critical behavior in the area difference will be observed. 
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In the higher dimension, various black objects admit the GL-like instabilities and their 
zero modes are expected to produce the new deformed branches. For example, Myers 
Perry black holes were conjectured to be unstable for the sufficiently large rotation [25], 
and it is numerically shown that each zero mode deformation leads to the new ‘bumpy’ 
solution [26, 27]. AdSs x S 5 black holes also have the deformed blanches [28]. These 
deformed solutions will show the similar critical behavior. The large D effective theory 
approach will be a good analytical tool to observe such nonlinear phenomenon. The authors 
have already shown that the bumpy Myers-Perry solution can be obtained by the 1/D 
expansion [19]. 

As the further extension of this work, it is interesting to study the time evolution of 
the non-uniformity in the 1/D expansion. This will be partially realized by including the 
slow time dependence, which is already formulated in [17]. For D > 13, the NUBS phase 
becomes stable and the time evolution of the GL instability can be followed by the current 
small gradient effective approach. Since the numerical simulation becomes more difficult 
in the higher dimension by the resolution problem, the large D analysis will be a suitable 
approach to this problem. For D < 13, the NUBS phase is no longer stable, and the GL 
instability develops to the self-similar structure which consists of the sequence of black holes 
and thin strings connected with each others [29-31]. Since the small gradient assumption 
will break down in such phase, we need to incorporate the large gradient (~ D/r o) dynamics 
which is not yet formulated. 

The large gradient dynamics will be also required to understand static solutions sup¬ 
ported by the interaction between the multiple horizons, because the small gradient degrees 
of freedom are confined in the near-horizon geometry. In section 6.3.1, we show that lo¬ 
calized black holes have the degenerate phase in the large D limit. This is because the 
gravitational force becomes non-perturvative in 1 /D, ( 2p$/L) D . To keep the deformation 
finite in the large D limit, one should assume very large black holes L /2 — po r*j Po/D, 
which will involve the large gradient dynamics. Near the merger point, the NUBS phase 
was conjectured to have the conical waist [32, 33]. In such phase, the gradient along the 
horizon will be also large. 
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A Uniform Black Strings in the conformal coordinate 

In this section, we will show the large D limit of the uniform black string in the conformal 
coordinate. D = n + 4 black string solutions are given by 

n n 

ds 2 = ~f(p)dt 2 + f(p)~ l dr 2 + p 2 dti 2 n+l + dZ 2 , f(p) = 1 - (A.l) 
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If the extra dimension is compactified with the circumference L, the mass and tension of 
the spacetime are given by 


M = 


(n + l)£l n+1 Lp% 


T = 


^n+lPo 


(A-2) 


167tG 167 tG 

where Q n +i = 27 r n / 2 + 1 /T(n/2 + 1) is the volume of S n+1 . This metric can be transformed 
to the conformal coordinate 

ds 2 = —F 2 (r)dt 2 + dr 2 + dZ 2 + p 2 (r)dQ^ l+l . 

by the transformation 

p -> r = p 2 Fi 

where r coincides with p at p —>• oo. F(r) is implicitly given by 

r 2-^1 


'1 

1 

i 

(»Y\ 

2 ’ 

1 

n 

n 

\p) _ 


r = 


Po 


-F{rf 

In n 


(1 - F(r) 2 )n 

In this coordinate, the horizon position F(rh) = 0 becomes 

(i -1) 


(A.3) 


(A.4) 


(A.5) 


rh = 


(A. 6 ) 


r(*-£) 

Large D limit Though eq. (A.5) cannot be solved explicitly, the large n limit allows the 
expression in the 1/n expansion. Keeping R = 4r n /pft finite, the limit n —> oo in eq. (A.5) 
leads to 

41 /n 


R i/n = 


(1 - F 2 )n 


T 2 


F\ 


111 , 

-,1-;1 - F 2 

2 ’ n n 


= 1 -tanh l F + 0(n ) 


n 


Comparing 0{n x ) terms in eq. (A.7), the leading term in F(R) is determined 

F(R) = ^ T \ + 0(n- 1 ). 

Recalling the relation F 2 ( R) = 1 — p 1 ^/p n , we also obtain 

i / 1 \ 2 /n 

p(R) = Po (1 - F 2 ( R))"» = roR 1 /" [l + -j + 0(n~ 2 ), 


(A.7) 


(A. 8 ) 


(A.9) 


where we introduced the horizon scale ro = 4 l ^ n po- Thus, in the conformal coordinate, 
the near-horizon geometry of the uniform black string in the leading order becomes 

2 / -, \ 4/n 

(A.10) 


ds 2 — — 


R- 1 


R + l 


dt 2 + dr 2 + dZ 2 + rg R 2/n ^1 + Fj 


dn 2 n+1 . 


We note that the horizon radius rh in eq. (A. 6 ) slightly differs to rg in the higher order, 


(-, C(2) 2C(3)\ 

rh-r o 1 -o-W ' 

\ n z r\r ) 


(A.ll) 
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B Near-zone solution 


B.l Next-to-leading order 

In the following, we use x = R/M(z) for the abbreviation. Li^(x) denotes the polyloga¬ 
rithm. 


= - 


4 M"(z)x(x — l)lnx 
M{z){x + l) 3 


+ 


4 (M(z) 2 — M'(z) 2 + M(z)M"(z))x(x — 1)(C(2) + In 2 x + 2 Li 2 (—x)) 

M(z) 2 (x + l) 3 


(B.l) 


<7« = 81n 2 (x 4 - 1 ) - i(M ' (Z) ~ <y M " {l)) + 1 ) 

M(z) 


+ ((3 + 5 x)M(z) 2 - (x - 1 )M'{z) 2 + (x - 1 )M(z)M"(z)) 


2 In 2 x 


(x + l)M(z) 2 


1 v J <- ■-i)m(z) y 

(C(2)x + (x - l)Li 2 (-x)), 


+ 


V ' ' ' (x + i)m( 2 ) 7 

4(M(z) 2 - M'{z) 2 + M{z)M"(z)) 


(x + 1 )M(z) 2 


(B.2) 


and 


+ 


(1) _ 8(M'(z) 2 — M(z)M"(z)) 2 2 

M(^)4 m + ij 

, 2 (M(z)M"(z) - 2 M’(z) 2 )(M(z) 2 - M\zf + M(*)M"(z)) „ 3 _ T . , „ 

+ 3 M(z )4 (in X + 5 

4(M(z)Af"(z) - M'{z) 2 )(M{z) 2 - M\zf + M{z)M"(z)) 
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( (1 — x\ , 2 C(2)x — 2x Li 2 (— x) 
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+ M(zf ( xM"(zf - (x + l)M'(z) 2 (21nM(z) - 1))). 
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B.2 Next-to-next-to-leading order 

The NNLO solution is given in the Mathematica filed attached. We only note that 
includes the following integration which we could not express by known functions 


Fi(x) = 


In 2 1 ln(t + 1 ) — In 3 1 
t(t- 1 ) 


dt = 


l 


x / x In 2 u\n{u + 1 ) 


1 — u 


du. 


(B.4) 


At x = 1, this integration converges to F\ (1) = 0.3457. Fortunately, we do not need this 
value for the current order analysis. 

C Next-to-Next-to-Leading order matching 

We show the matching condition up to 0(n~ 2 ). The detail of Aq(z) is given by 
4C(2)M(*) - 4C(2 ^) (g)2 + 4(C(2) - In M(z))M"(z) 


Ao(z) = 4 M(z) + 
2 


n 


+ 5 n 2 M(z) 3 [ 5 M (A > 4 (C(2 ) 2 + 4C(3)) - (-30C(2) + C(2 ) 2 + 20 In M(z) - 5((3)) M'(z) 4 

+ M(z) 3 (20C(2) + 7C(2 ) 2 + 10(-1 + 2C(2))ln M(z) - 15 In 2 M(z) + 25C(3)) M"(z) 

- M(z ) (50C(2) + C(2 ) 2 - 40 In M(z) + 51n 2 M(z) + 10C(3)) M'(z) 2 M"(z) 

- M(z) 2 ((10C(2) + 4C(2 ) 2 + 20C(2) In M(z) + 10 In 2 M(z) + 25C(3)) M\z ) 2 

- (20C(2) + 2 C( 2) 2 - 10 In M(z) + 10 In 2 M(z) + 5C(3)) M"(z) 2 )]. (C.l) 
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With eq. (C.l) and Bq(z), eq. (4.22) becomes 
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where the higher derivatives of M(z) are already eliminated by using eq. (4.20). The 
nonlinear terms of M"(z) in the higher order can be also eliminated by the repeated use of 
the same equation. After the integration and again eliminating M'(z ) in the higher order, 
we obtain eq. (5.1). 


D Near uniform solution 

The solution of eq. (5.1) expanded with 5=1 — e 2 /2 is expressed in the Fourier series 

OO 

M(*) = 5^J“j( e ) eJ cos (jk(e)z) (D.l) 

j =0 
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where the wavenumber k(e) is given by eq. (5.20). In the following, we show first few terms. 
The terms up to 0(e‘) are given in the attached Mathematica file. 
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